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ABSTRACT 





The non-availability of ellipsoidal heights of local geodetic Datums has made it necessary for the 
application of ellipsoidal heights transformation models to the available global ellipsoidal heights to 
obtain their respective theoretical heights in local Datums. It is required to know the accuracy, as 
well as reliability of any model of interest before its application. For that reason, this study 
comparatively analyses the Molodensky and Kotsakis models for the transformation of ellipsoidal 
heights between geocentric and non-geocentric Datums to determine the reliability of the Kotsakis 
model. The Global Navigation Satellite System (GNSS) data of the used stations were processed in 
World Geodetic System 1984 (WGS84) datum to obtain their global geographic coordinates and 
ellipsoidal heights. The coordinates, ellipsoidal heights and the transformation parameters between 
WGS84 and Minna Datums were applied to the Molodensky and Kotsakis models to compute the 
Clarke 1880 theoretical heights of the stations. The Molodensky model was used as a reference to 
which the Kotsakis model ellipsoidal heights were compared to obtain the Kotsakis model 
ellipsoidal heights discrepancies, as well as residuals. The residuals were used to compute the 
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Root Mean Square Error (RMSE) of the Kotsakis model. The computed RMSE, as well as reliability 
of the model is 1.244 m. The study concluded that the low reliability, as well as accuracy of the 
Kotsakis model might be as a result of the two rotation datum shift parameters in it as they are the 


main differences between the two models. 
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1. INTRODUCTION 


Practical height computation from the processed 
observed Global Navigation Satellite System 
(GNSS) data requires the application of geoid 
height to the ellipsoidal height obtained from the 
processed observations. Ellipsoidal heights are 
theoretical heights obtained from the GNSS 
observations, which are measured from the 
surface of the reference ellipsoid to the observed 
point on the earth surface [1]. Non-geocentric 
datum ellipsoidal heights are not readily available 
in most of the GNSS observation areas and 
regions. Most GNSS height adjustments, as well 
as fitting during observations processing are 
done using the orthometric heights of the existing 
controls, used as reference stations in the 
observation. However, ellipsoidal heights are 
applied to the GNSS observation for theoretical 
height adjustment, likewise orthometric height to 
spirit levelling for practical heights reduction. The 
orthometric heights are measured along the 
gravity vector direction and referenced to the 
geoid, as well as the mean seal level [2]. The 
ellipsoidal and orthometric heights have their 
respective reference surfaces. The erroneous 
use of the orthometric height for the GNSS 
observations processing to obtain local 
ellipsoidal heights of points is as a result of the 
unavailability of ellipsoidal heights in the 
observation area or region. The Clarke 1880 
ellipsoid adopted for geodetic computation in 
Nigeria is flatter and bigger compared to the 
WGS 84 spheroid (see Fig. 1). Nigeria is located 
between latitudes 4°N and 14°N, which is closer 
to the equator than the North Pole. So, it is 
expected that the ellipsoidal heights computed 
on the Clarke 1880 are smaller in value than 
those computed on the World Geodetic System 
1984 (WGS84) ellipsoid. 


The local, as well as non-geocentric datum 
ellipsoidal heights can be obtained from the 
conversion of ellipsoidal heights computed on the 
global, as well as geocentric (WGS84) ellipsoid. 
The conversion can be achieved through the 


application of the 5-parameters Molodensky’s 
model [8,4,5] and the 8-parameters Kotsakis 
model [6] for ellipsoidal heights transformation 
between the geocentric and non-geocentric 
Datums, as well as reference frames. The 
Molodensky model involves the use of the 3 
translation datum shift parameters, change in 
semi-major axis and difference in flattening 
between the two reference frames, as well as 
ellipsoids while the Kotsakis model comprises 3 
translation and 2 rotation datum shift parameters, 
change in scale, change in semi-major axis and 
difference in flattening between the two reference 
ellipsoids. The Molodensky method was recently 
used by [7] for the determination of the ellipsoidal 
height of the Nigerian geodetic/Minna datum and 
Root Mean Square Error (RMSE) of 0.00m was 
achieved. The method was compared with other 
two methods and was recommended as the best 
among the three methods. The Kotsakis method 
has not really been applied to Nigeria. Here, the 
Molodensky method is used as a reference to 
which the Kotsakis method is compared to 
determine its reliability. Consequently, this study 
comparatively analyses the Molodensky and 
Kotsakis ellipsoidal heights transformation 
between geocentric and non-geocentric Datums 
models to determine the reliability of the Kotsakis 
model. 


2. METHODOLOGY 


The adopted methodology involves the 
transformation of geocentric datum (WGS84) 
ellipsoidal heights obtained from the GNSS 
observations to local ellipsoidal height in the 
Nigeria Minna datum using the Molodensky and 
Kotsakis methods and comparing their results. 
The application of the two methods requires the 
use of the 5-parameters Molodensky change in 
ellipsoidal height computation model, 8- 
parameters Kotsakis model, datum shift, as well 
as transformation parameters between the 
WGS84 and Minna Datums, and the two 
Datums, as well as ellipsoids properties (semi- 
major axis and flattening). 
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Fig. 1. Relationship between WGS84 and Clarke 1880 ellipsoids 


2.1 The 5-Parameters Molodensky Model 


The 5-parameters Molodensky model used for the transformation of ellipsoidal heights between 
geocentric and non-geocentric reference frames is [3,4,5,7,8]. 
b , 
J>a(2)e. sin’ p 
a 


a 


Ahyosg4 = Ty cos pcos A+T, cosgsindA+T, sino bf (1) 


N 


Where, 


T,,, 1, T; =Translation parameters between WGS84 and Minna Datum. 


gy, A = Geographic coordinates (Latitude and Longitude) of points. 


a = Equatorial radius of the Clarke 1880 ellipsoid. 
b= Polar radius of the Clarke 1880 ellipsoid. 
f = Flattening of the Clarke 1880 ellipsoid. 


Aa = Change in equatorial radius between the two ellipsoids (Minna minus WGS84) 
Af = Change in flattening between the two ellipsoids (Minna minus WGS84) 


a 
(1-e? sin? oy? 


2 2 
a’ —b 
2 
a 


R, = Radius of curvature in prime vertical = 








e= Eccentricity squared = 2 f —f? oa 


Nctarkei 880 7 Nyassa + Alwesga 


b=a(l-f) (4) 


2.2 The 8-Parameters Kotsakis Model 


Having computed Afy,.g,, the non-geocentric 
(Clarke 1880) ellipsoidal heights (/ 


Clarke 880) are 
obtained using [5,7] 


h'-h=oh(t,.) + On, ) + ONL) + ONCE, 


+ 0h(da) + dh( oF ) 
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The 8-parameters Kotsakis model used for the 
transformation of ellipsoidal heights between 
geocentric and non-geocentric reference frames 
is [6]. 
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Where, 
ON(t,.) =t,, cos pcos A (7) 
ONL, ) =t, cos psin A (8) 
ONn(t,) =t, sin@ (9) 


OMé,) =—é,Ne* singcos gsin A 
n(é,) = €,Ne’ sin pcos pcos A 
on( ds) = (aW + h)ds 


0h(da) = —Woéa 


Wf) = AO sin? pf 


W=,l-e’ sin’ 9 


N in equations (10) and (11) is the radius of 
curvature in prime vertical as given in equation 


(2) 


The quantities da=Aa=a'-—a and 
of = Af = f'—f correspond to the difference 
in the numerical values for the semi-major axis 
and the flattening of the reference ellipsoid, as 


these are used in the respective reference 
frames, GRF1 and GRF2 [6]. 


(10) 


2.3 Transformation Parameters between 
WGS84 and Minna Datum (Clarke 
1880 Ellipsoid) 


The transformation parameters from WGS84 to 
Minna datum are [9] 
T, = 93.809786m + 0.375857310m 
T, =89.748672m +0.375857310m 
T, = —118.83766m + 0.375857310m 
Ry, =0.000010827829 + 0.00000103 11322 
R, = 0.0000018504213 + 0.0000015709539 
R, = 0.0000021194542 + 0.0000013005997 
S = 0.99999393 + 0.0000010048219 

hh 


2.4 The Nigerian Geodetic and WGS 84 
Datums 


The Nigerian geodetic datum (Clarke 1880 
ellipsoid) and WGS84 ellipsoid semi-major axes 
(a) and flattening (f) are respectively [10] 
6378249.145 m and 1/293.465, and 
6378137.000 m and 1/298.257223563. 


da = Aa =a'—a=112.145 
ff =Af = f'— f =0.0000547507139518535 
2.5 Root Mean Square Error (RMSE) 


The Root Mean Square Error (RMSE) of a model 
is computed to indicate its accuracy, as well as 
reliability. Here, the RMSE is computed by 
comparing the transformed ellipsoidal heights 
obtained from the two models using the 
Molodensky model as aé_ reference. The 
computation of the RMSE of the transformation 
model is done using [11]. 


RMSE = LS 6? 
N jt 


Where, 


(12) 


2 _ 2 
oS (ypotodensky ae gotsaris) 


= Molodensky model ellipsoidal height 


Nyotodensky _ 


Kotsakis = KOtsakis model ellipsoidal height. 


n= Number of points. 


A total of 11 GNSS points located within Edo 
State were used in the study. The observation of 
the points was carried out with 5 dual- 
frequencies GNSS receivers. The geographic 
coordinates and ellipsoidal heights of the points 
were processed on the WGS84 ellipsoid using 
the Compass Post-processing software as the 
study involves the transformation of global 
dataset to local. Table 1 shows the geocentric 
(WGS84) datum geographic coordinates and 
ellipsoidal heights of the used stations. 


The changes in ellipsoidal heights between the 
WGS84 and Clarke 1880 spheroids and the 
Clarke 1880 ellipsoidal heights regarding the 
Molodensky model were respectively computed 
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using equations (1) and (5) while those of the Figs. 2 and 3 for Molodensky and Kotsakis 
Kotsakis model were computed using equation models respectively. The reliability, as well as the 
(6). The computations were done with computer root mean square error of the Kotsakis model, 
programs developed in this study, as shown in was computed using equation (12). 


















































































































































































































































Home Insert Page Layout Formulas Data Review View PDFelement 
: = * 
& cut Calibri “1 Bl & | Fwep tet General . Ee 7 i = aK = AutoSum Ay A 
Ba Copy + Fill + 
Paste BsIuU- z necaess Comen + % 9 | 32 399 Conditional Format Cell | Insert Delete Format Sort & Find & 
+ GF Format Painter ia = Eid Mera $ a “0° © | Formatting + as Table Styles ~ . . . @ Clear > Filter Select 
Clipboard ey Font Fy Alignment Ey Number Fy Styles Cells Editing 
139 qc fe | =-118.83766*(SIN(D39*0.0174532925199432)) “ 
= 
eer B ce D E F G ay J K L M N ° P Cl 
1 
,. EXCEL PROGRAM FOR TRANSFORMATION OF GLOBAL ELLIPSOIDAL HEIGHTS TO LOCAL (DECIM. 
seoos“Eccentricity 0.006803511 = _ 6378137,000 = 0,003352811 _ == 6378249.145 == 0,003407561 
3 biome y awess4 fwoss4 a f Minna 
4 
5 2 F Aa ad 
STATION | LATITUDE LONGITUDG ° a Neoew | [2V@.) aNt,) Nt.) e Ww Nev a = or, sng?) | |Afigesss bya 
5 x = 
6 PBG134| 6.649502 [6.4525807| 6.649591731 | 6452589697] 288.3613 | 92.58844048 | 10.0181983 | -13.76103384 | 0006694381 | 0999955117 | 6378423.284 | -112.139966569796| 4666979063972 | _-18.6274 269.7339 
7 _PBG135 | 6.649729 [6.4536393| _6.649728033 | 6.453639294| 287.9072 | 92.58823179 | 10.0198182 | -13.7613165 | 0.006694381 | 0999955115 | 6378423.296 | -112.139966363012| 4667170797990 | _-18.6261 269.2811 
8 PBG137 | 6.668297 |6.5692762| 6.668296786 | 6.569276181| 79.7948 | 92.56340358 | 1019818791 | -13.79956843 | 0.006694381 | 0.99995a865 | 6378424.800 | -112.139938339778| 4693154432278 | _-18.4848 61.3100 
° PaG138 | 6.66879 | 6.569233 | 6.668789994 | 6.569233031| 79.326 | 92.56331844 | 10.19811095 | -13.g00584a8 | 0.006694381 | 0.99905ae58 | 6378424932 | -112.139937594366| 4693845591842 | -18.4852 60.8408 
10 PBG139 | 6.668562 | 6569038 | 6.668562142 _| 6.569937992| 75.3667 _| 92.56323032 | 10.19920528 | -13.80011509 | 0.006694381 | 0.999954861 | 6378424913 | -112.139937938738| 4693526284033 | -18.4841 56.8826 
11 BEMB06 | 6.301587|5.6311678| 6301587494 | 5.631167753| 97.4262 _| 92.79300446 | 8753313072 | -13.0438415 | 0.006694381 | 0.999959673 | 6378394.219 | -112.140477568832| 4193172930335 | -19.4448 77.9814 
a2 ENT" | 6.302732 |5.6311172| _6.302732483__| 5.631117192| 96.5354 _| 92.79280774 | 6.753215413 | -13.04620198 | 0.006694381 | 0.999959659 | 6378394312 | -112.140475931847| 4194690764632 | _-19.4460 77.0894 
13 EDRPO1|6.078587|5.6680006| 6.078586636 | 5.668000589| 58.581 | 92.82628065 | 8.814i01198 | -12.58401053 | 0.006694381 | 0.999962467 | 6378376402 | -112.140790810457| 3902731564757 | _-19.1817 39.3993 
14 EDRPO2| 6.078634 |5.6696422| 6.078633864 | 5.669642153| 58.5242 | 92.82600851 | 8.816644825 | -12.58410794 | 0.006694381 | 0.999962466 | 6378376406 | -112.140790745294| 3902791984644 | _-19.1795 39.3447 
15 UHA100 | 6.740426 [64317531] _6.740425822 _| 6.431753056| 188.3327 | 9257502204 | 9.984125673 | -13.9ae14ae4 | 0.006694381 | 0.999953ee8 | 6378431123 | -112.139828752342| 4794765888309 | _-18.7341 169.5985 
6 _UHA101 | 6.740359 |6.4610593| 6740359064 _| 6.461059308| 184.7541 | 9256968481 | 10.02942712 | -13.94801134 | o.006694381 | 0999953889 | 6378431117 | -112.139828854304| 4794671347069 | _-18.6941 166.0600 
7,0 0 oO 0 o 0 93.809786 o 0 (0.006694381 1 6378137.000_| -112.144999999552| 0.000000000000 | __-18.3352 -18.3352 
is, 0 o 0 oO 0 93.809786 ° 0 0,006692381 1 6378137,000_| -112.144999999552| 0.000000000000 | _-18.3352 -18.3352 
wl 0 O oO 0 0 0 93.809786 0 0 10.006694381 1 6378137.000_| -112.144999999552| 0.000000000000 | _-18.3352 718.3352 
2,0 O 0 0 0 0 93.809786 0 ° (0.006694382 1 6378137.000_| -112.144999999552| 0.000000000000 | __-18.3352 "18.3352 
ul 0 0 0 0 0 O 93.809786 O ° (0.006694381 1 6378137.000 | -112.144999999552| 0.000000000000 | _-18.3352 "18.3352 
210 0 i ers a = 93,809786 ry ry 0,008694381_| 1 6378137.000_|-112.1449999995521 o.o00000000000 | -18.3352 = ud 
i «> || PROCESSING TAB 4TERPOMEEM Sheets el Dal i 0 





Fig. 2. Molodensky model change in ellipsoidal heights computation 
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Fig. 3. Kotsakis model Clarke 1880 ellipsoidal heights computation 
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3. RESULTS AND DISCUSSION 


Table 2 presents the discrepancies in the 
ellipsoidal heights and Root Mean Square Error, 
RMSE of the Kotsakis model. They were 
computed to show the range of _ the 
discrepancies, as well as the differences in 
ellipsoidal heights between the two models and 
the accuracy of the Kotsakis model relative to the 
Molodensky model. It can be seen from Table 2 
that the minimum and maximum discrepancies of 
the Kotsakis model ellipsoidal heights are - 


0.1775m and 1.7459m respectively. It implies 
that the Kotsakis model ellipsoidal heights 
discrepancies range from -0.1775 to 1.7459m. 
The computed range is limited to the used 
stations, as well as the location of the points 
(Edo State). | can also be seen in Table 2 that 
the RMSE of the Kotsakis model is 1.244m which 
implies that the model has a reliability, as well as 
accuracy of 1.244m. The low accuracy of the 
model might be as a result of the two rotation 
datum shift parameters terms in it since they are 
the main differences between the two models. 


Table 1. Geographic latitudes, longitudes and ellipsoidal heights of stations 











Station WGS 84 Datum 

Latitude (Decimal. Longitude (Decimal. Ellipsoidal 

Degree) Degree) Height (m) 
PBG134 6.649591 731 6.452589697 288.3613 
PBG135 6.649728933 6.453639294 287.9072 
PBG137 6.668296786 6.569276181 79.7948 
PBG138 6.668789994 6.569233031 79.3260 
PBG139 6.668562142 6.569937992 75.3667 
BEM606 6.301587494 5.631167753 97.4262 
ENV100D 6.302732483 5.631117192 96.5354 
EDRP01 6.078586636 5.668000589 58.5810 
EDRP02 6.078633864 5.669642153 58.5242 
UHA100 6.740425822 6.431753056 188.3327 
UHA101 6.740359064 6.461059308 184.7541 





Table 2. Transformed ellipsoidal heights discrepancies and RMSE of Kotsakis model 











Station NeLoBa Atocav/minna (M) Difference/ Difference 
Molodensky Kotsakis Discrepancy (m) Squared (m?) 
Model Model 
PBG134 288.3613 269.7339 268.2733 1.4606 2.1335 
PBG135 287.9072 269.2811 267.8179 1.4633 2.1411 
PBG137 79.7948 61.3100 59.5655 1.7446 3.0436 
PBG138 79.326 60.8408 59.0971 1.7436 3.0402 
PBG139 75.3667 56.8826 55.1367 1.7459 3.0482 
BEM606 97.4262 77.9814 78.1566 -0.1752 0.0307 
ENV100D 96.5354 77.0894 77.2670 -0.1775 0.0315 
EDRP01 58.581 39.3993 39.0484 0.3509 0.1231 
EDRP02 58.5242 39.3447 38.9894 0.3554 0.1263 
UHA100 188.3327 169.5986 168.3520 1.2466 1.5540 
UHA101 184.7541 166.0600 164.7334 1.3266 1.7599 
Kotsakis Model RMSE = 1.2443m 





4. CONCLUSION 


The study has comparatively analyzed the 
Molodensky and Kotsakis ellipsoidal heights 
transformation between geocentric and non- 
geocentric Datums models and determined the 
accuracy, as well as reliability of the Kotsakis 
model. The study has determined the range of 
the discrepancies of the Kotsakis model limited 
to the used stations to be -0.1775 to 1.7459m. It 


has also determined the accuracy of the Kotsakis 
model to be 1.244m. It again stated that the low 
accuracy of the model might result from the two 
rotation datum shift parameters terms in it. 
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